Specific heat anomaly in the d-density wave state and emergence of incommensurate 

orbital antiferromagnetic order 
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We study the effect of finite chemical potential on the d-density wave state that has been proposed 
to explain the pseudogap phenomena in underdoped cuprates. We find that the specific heat anomaly 
at the transition temperature, below which the d-density wave state forms, gets weaker when finite 
chemical potential is introduced. This provides a useful ground for the proper interpretation of 
the specific heat measurement in regard to the existence of the d-density wave state below the 
pseudogap temperature. Further increase of the chemical potential leads to an incommensurate 
orbital antiferromagnetic state before the system eventually turns into the normal state. This is 
an inhomogeneous state characterized by novel charge ordering and an analog of the Fulde-Ferrel- 
Larkin-Ovchinnikov (FFLO) state in superconductivity. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.10.Hf 



a. Introduction One of the recent interesting pro- 
posals made in high temperature superconducting 
cuprates is that the so-called d-density wave state is 
responsible for the pseudogap phenomena. || The d- 
density wave state is a particle-hole paired state with a 
finite momentum Qo — (tTjTt). The particle-hole pair 
wavefunction has a d-wave symmetry, /(k) cx cosfc^ — 
cos ky , that occurs due to the formation of staggered 
currents circulating in the square lattice. Therefore, the 
d-density wave state breaks the translational and time 
reversal symmetries. ||, ^ 

Since the d-density wave state is an orbital antiferro- 
magnetic state, there should be an elastic neutron scat- 
tering signal at the momentum Qp. This is a direct con- 
sequence of the broken symmetry. It was argued that the 
recent neutron scattering experiment is consistent with 
the existence of the d-density wave state in underdoped 
cuprates ^. On the other hand, there is no clear 
experimental evidence of the broken symmetry state in 
thermodynamic quantities, such as specific heat at T* 
where the pseudogap develops. It has been regarded that 
this is at odds with the proposal, because one may expect 
that there should be a clear specific heat jump due to the 
translational symmetry breaking. In fact, in the case of 
the half-filled band, the normalized specific heat jump for 
the d-density wave state is the same as that of the d-wavc 
superconducting state within the mean field theory. This 
is simply because these states have the same form of the 
quasiparticle spectrum, E]^ = + Aj^ when the chem- 
ical potential is zero {fi — 0). Here Ak = Aq cos (20k) 
is the d-wave gap in the d-density wave state or in the 
d-wave superconductor. The specific heat measurements 
show a clear evidence of the phase transition at Tc, 
while there is no anomaly reported at T*. 

This conclusion should be, however, reexamined when 
the system is away from half- filling. In the presence of 
the finite chemical potential, /i, the quasiparticle spec- 
trum is quite different for two cases. In the d-wave 
superconductor, the quasiparticle dispersion is given by 
i^k = \/ (ck — m)^ + A^, while in the d-density wave state 



it is given by i?k = — At± V^k+^k- The difference comes 
from the fact that one is a particle-particle paired state 
and the other is a particle-hole paired state. 

In this paper, motivated by the above observations, we 
study the effect of the finite chemical potential on the 
commensurate orbital antiferromagnetic state dubbed 
"d-density wave state". Using a BCS-like Hamiltonian 
for the d-density wave state, we show that 1) the specific 
heat anomaly at the transition temperature, T = T^, due 
to the d-density wave ordering gets weaker as the chemi- 
cal potential increases. Therefore, when the chemical po- 
tential is appreciably large, it will be hard to determine 
the existence of the d-density wave state by measuring 
the specific heat. In view of the importance of the fi- 
nite chemical potential, we ask whether there is a phase 
transition driven by the chemical potential. We find that 
2) the stability of the d-density wave state gets weaker 
as the chemical potential increases, and beyond a critical 
value, an incommensurate orbital antiferromagnetic state 
occurs in a range of chemical potential. It is an inhomoge- 
neous state characterized by charge ordering in contrast 
to the d-density wave state with uniform charge density 
and an analog of the FFLO state in superconductivity. ||] 
Certain choices of the incommensuration correspond to 
the stripe patterns of the charge ordering similar to those 
discussed in the context of the spin density wave in the 
cuprates. Q 

b. Role of the chemical potential The d-density wave 
state is a particle-hole paired state where the band dis- 
persion of the particles and holes satisfy ek — — ek-i-Qo 
with Qo = (vr, n). Within mean field approximation, the 
d-density wave state can be stable in an extended Hub- 
bard model at half-filling, i.e., = 0, where the particle- 
hole symmetry exists. Q 

When the chemical potential becomes finite, /i ^ 0, 
the number of particles is not equal to that of holes. 
Therefore, the chemical potential suppresses the ten- 
dency of paring and effectively acts as a pair breaker. 
In other words, the nesting condition gets weaker with 
finite chemical potential, so does the stability of the d- 
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density wave state. In fact, the chemical potential can 
be regarded as an analog of the Zeeman magnetic field in 
superconducting states, where the Zeeman field splits the 
up- and down-spin bands and acts as a pair breaker. We 
will show below that the gap equation of the d-density 
wave state with finite /i, Eq.|^, is the same as that of 
the d-wave superconducting state with a finite Zeeman 
field. [|lO[ Q Once we realize this simple analogy, it is 
clear that the specific heat anomaly will be reduced by 
introducing the chemical potential. This is because the 
Fermi surface is no longer a point, but there will be a hole 
pocket whose size depends on the magnitude of fi. There 
were similar studies in the spin density wave showing that 
the specific heat jump gets smaller when the nesting con- 
dition gets weaker. [|l2| 

c. Specific heat anomaly in the d-density wave state 
In order to compute the specific heat jump at the tran- 
sition temperature, T^, we use the following BCS-like 
Hamiltonian. 

k,cr 

+ 5kk'Q cl ,^Ck+Q,aCk'-QM'cl, ,^,, (1) 

k,k',Q,cT,a"' 

where Ck is the band dispersion and satisfies ek = 
^Ck+Qo- 5kk'Q is the pairing potential that takes the 
following form when Q — Qq. 

5kk'Qo = -5cos(26'k)cos(26ik')- (2) 

The gap equation is obtained in the mean field approx- 
imation as follows. 
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where 



e^-t-A^ cos^ 26* and we define /i to be positive. 
Here A^(0) is the density of state at e = 0, and we have 
assumed that there is no singularity in A^(0).||l^ 

Using the standard procedure, one can obtain the fol- 
lowing formula for the transition temperature Td{fJ.) at 
finite chemical potential. 

ln^=RcV.fi + ^V^fiV (4) 
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where T^o is the transition temperature when fj, = and 
ip{x) is the di-gamma function. As expected, Td{^) de- 
creases as the chemical potential increases. 

Now let us evaluate the specific heat jump at T — Td- 
Since we are interested in the specific heat anomaly at 
T = Td, we expand the gap equation in A(u, T) / ttT. Us- 
ing the conventional method described in [Q , the specific 
heat jump at Td can be computed as 



CddM) - CN{Td) = N{0) 
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where C,{s,a) = I]^^o ^/('^ + ^^sre Cddw and Cn 
represent specific heat of the d-density wave state and 
normal state respectively. 

The specific heat coefficient ^ddw — Cddw IT cit T — Td 
is shown in Fig. 1 as a function of fJ-/Td, where in the 
normal state is set to be unity {'^n — 1). As shown in Fig. 
1, the specific heat anomaly, d^ddw/lN = {,lddw-lN)llN 
gets smaller as /i increases. Now let us compare the ab- 
solute value of the specific heat anomaly, S^ddw/^N, with 
that of the d-wave superconducting state, S^sc/lN ~ 0.9. 
In the case of the d-density wave state, when fi — O.GT^ = 
0.56Tdo, 5-fddw/lN = 0.26. 
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FIG. 1: The specific heat coefficient, 'yddw/lN is shown as a 
function of ^/Td- The dashed line represents 7jv in the normal 
state, set to be unity. Notice that ^i/Td is a monotonically 
increasing function of ^/Tdo- 

Thus, the specific heat jump associated with the d- 
density wave state can be quite small compared with that 
of the d-wave superconducting state when the chemical 
potential is finite. Notice that our results are obtained 
in the mean field theory and the anomaly can be even 
smaller once we take into account the fiuctuation effects. 

d. Phase transition Provided that the chemical po- 
tential in the d-density wave state plays the same role as 
the Zeeman magnetic field in superconductors, let us first 
review the phase transition driven by a Zeeman field in 
superconductors. It is clear that sufficiently strong Zee- 
man field will eventually drive a phase transition into 
the normal state. The nature of the transition is de- 
termined by the competition between two energy scales; 
the condensation energy ^N{0)A^ and the magnetic en- 
ergy where x is the spin susceptibility proportional 
to the density of state N{Q). This competition leads to 
the phase transition between the BCS state and normal 
state at a critical Zeeman field He = A/(\/2/iB) (/is is 
the Bohr magneton) for s-wave superconductors. On the 
other hand, it has been also known that, for a range of 
intermediate strength of the Zeeman field around He, it is 
more favorable to pair up- and down-spin electrons across 
the spin-split Fermi surfaces and form an inhomogeneous 
superconducting state. This is called the Fulde-Ferrel- 
Larkin-Ovchinnikov (FFLO) state. Il^, The phase 
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transition from the BCS superconducting state to the 
FFLO state is the first order transition for both s-wave 
and d-wave superconductors |ll|, [Tst - 

The d-density wave state with finite chemical poten- 
tial should experience the same kind of phase transition 
due to the existence of two competing energy scales. The 
condensation energy favors the particle-hole paired state. 
As the chemical potential increases, however, the normal 
state is favored by iV(0)/i^. Therefore, the qualitative 
phase diagram should be similar to those obtained in 
superconductivity pH, pT[ . What is really interesting is 
that there is a window of chemical potential where a novel 
phase, an analog of the FFLO state, is favorable over 
the commensurate orbital antiferromagnetic and normal 
states. This phase is the incommensurate orbital antifer- 
romagnetic state where the particle-hole pairing occurs 
with a momentum Q — Qo + q, where \q\ oc fi at T — 0. 
Similar studies were undertaken in the context of the ex- 
citonic insulator and an inhomogeneous phase with the 
s-wave symmetry was identified, p^ ] 

e. Incommensurate orbital antiferromagnetic state 
The incommensurate orbital antiferromagnetic order 
arises when a staggered current is circulating in the plane 
but its periodicity is incommensurate with the underly- 
ing square lattice. 



->- 



1.5/\ ^ \/-^ /\ \/ 



->- 



->- 
1 

->- 



->- 



-1.5\/-^ /\ ^ \/ /\ 



-1 



->- 



->- 



/N t /N t 



(a) 



9-o-m-o 

0- 9-0-9 

1 I I I 

(b) 



FIG. 2: (a) An example of the incommensurate orbital an- 
tiferromagnetic state with q = (0,7r/3). The numbers repre- 
sent the relative strength of each bond current. Since Qx ~ 0, 
the current is alternating with the same strength along the 
x-direction. (b) The charge modulation in real space for 
q = (0,7r/3). The black and white circles represent the am- 
plitudes of the charge modulation proportional to sin(7r/3) 
and — sin(7r/3) respectively. There is no charge modulation 
at the empty sites. 



Current modulation: It is straightforward to ob- 
tain the expectation value of the bond current in the in- 
commensurate orbital antiferromagnetic state, (Jrr') — 



current modulation is obtained as 

{Jr,r+x - Jr,r+y) » y^(CkCk+Qo+q) 
k 

X { [/(k, q) - ig{K q)] cos [(Qo - 
-[g(k,q) + i/(k,q)]sin[(Qo4 

Here /(k, q) and g(k, q) are given by 
/(k, q) = cos fc:r — cos /Cj, 
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+ cos {kx + qx) - cos {ky + qy), 
g(k, q) = sin fcj, — sin /ca; 

+ sin {kx + qx) - sin {ky + qy). 



(7) 



The following order parameter of a particle-hole paired 
state with the momentum Q = Qo + q is the simplest 
choice for the incommensurate orbital antiferromagnetic 
order. 



CkCk+Qo+q> 



4^Qo+?/(k,q), 



(8) 



cJ,Cr). Taking r' = y + x and r' = r + the 



where /(k, q) is given by Eq. (^. This order parameter 
of the particle-hole pairing is an analog of the FFLO 
order parameter of the particle-particle pairing. It is also 
reduced to the order parameter of the d-density wave 
state when Q — Qo. 

The order parameter, Eq. ^, leads to the following 
current pattern in real space. 

(Jr,r') = (-1)'"Aq„+, 

X [5x' ,x+i5y> ,y{cos (q • r) + cos (q • r + q^)} 

+ 5x',x-i5y',y{cos (q • r) -I- cos (q • r - q^)} 

- 5y,^y+i5x',x{co'i (q ■ r) + cos (q • r + qy)} 

- 5y' ,y-i5x' A^os (q • r) + cos (q • r - qy)}]{Q) 

One example of the incommensurate staggered current 
patterns is shown in Fig. 2 (a), where q = (0, tt/S). 
Different choices of the ordering vector Q lead to different 
patterns of the incommensurate staggered currents. For 
q — (0,7r/3), the current density has the modulation 
with the period [2a, 3a] in the x- and y-directions, where 
a is the lattice constant. Notice that there exists a set of 
the vanishing bond currents, Jr,r+y = 0, that are aligned 
along the a;-direction. This "bond-centered" stripe-like 
structure leads to the charge modulation (computed in 
the next section) that has the same periodicity as that of 
the current, as shown in Fig. 2 (b). 

On the other hand, when q — (0, 7r/4), a set of the 
vanishing bond currents in the a;-direction, Jr,r+x = 0, 
exists. This corresponds to the "site-centered" stripe-like 
structure along the a;-direction in the current pattern. As 
a result, the periodicity of the current is [2a, 8a], while 
that of the charge density is [2a, 4a]. That is, the period- 
icity of the current along the y-direction is twice larger 
than that of the charge. 

Charge ordering in the incommensurate orbital 
antiferromagnetic state: The staggered current pat- 
tern in the d-density wave state gives no charge modu- 
lation due to X]kGi?Bz(cos fcx — cosfcj^) = 0, where RBZ 



4 



stands for the reduced Brillioun zone. This property dis- 
tinguishes the d-density wave state from the ordinary 
charge density wave state. 

In the incommensurate orbital antiferromagnetic state, 
we find that there exists a charge modulation. As an 
example, let us compute the charge density when the 
staggered current pattern is given by Fig. 2 (a), where 
q = (0, 7r/3). The local charge density can be obtained 
from 

(^t(r)V-(r)) - ^^(4ck') exp [*(k - k') • r] 

k k' 

= no + AQ„+qh{qy) sin [(Qo + q) • r] , (10) 

where 

^ilv) — ^ [2cosfca; — cosfcj^ — cos (fcj^ + (7j^)] (11) 

The resulting charge modulation is shown in Fig. 2 (b). 
Notice that the stripe patterns of the charge ordering oc- 
cur for q = (0, ±|gy|) and q = (±|(72,|, 0). The amplitude 
of the charge modulation depends on the value of qy via 
h{qy) (for example, /i(7r/5) = —0.42) and it vanishes at 
qy = 0, which confirms the absence of a charge modula- 
tion in the d-density wave state. A similar, but different, 
relation between the charge density wave and the spin 
density wave was discussed in [| and ||. The main dif- 
ference between the stripe-like structure found here and 
those discussed in the context of the spin density wave is 
that there exists an "internal" charge modulation inside 
each charge stripe in the present case. 

A checker-board pattern of the charge ordering can 
also occur if q = (±|qj:|, ±|gy|). The checker-board pat- 



tern gets distorted as the asymmetry between q^ and qy 
increases. 

/. Discussion and Conclusion We show, at the 
mean-field level, that the specific heat jump at the transi- 
tion temperature of the d-density wave order gets weaker 
when the chemical potential is finite. As a result, it be- 
comes more difficult to determine the existence of the d- 
density wave state as the system is away from half-filling. 
The relevance of our results to the cuprates depends on 
the doping dependence of the chemical potential, which 
is currently not well understood. There is an indica- 
tion in La2_xSra;Cu04 that the chemical potential may 
be pinned in the underdoped regime. In this case, 
our result of the weak specific heat anomaly does not 
apply. Our results are more relevant when the chemi- 
cal potential moves sensitively as the doping concentra- 
tion changes, which seems to happen near the optimally 
doped regime. ||20|] Further increase of the chemical poten- 
tial leads to a phase transition into the incommensurate 
orbital antiferromagnetic state, where a novel charge or- 
dering occurs. The charge ordering patterns depend on 
the nature of the incommensuration of the current den- 
sity and can be observed by X-ray scattering. Some of 
the possible charge ordering patterns have stripe struc- 
tures similar to, but different from, those discussed in the 
context of the spin density wave in the cuprates. [|[ ^ 
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